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LONG-TERM  GOAL 

The  long-term  goal  of  this  project  is  to  improve  our  ability  to  understand,  model  and  predict  lateral 
mixing  and  the  associated  submesoscale  physical  structure  and  processes  in  the  upper  and  interior 
ocean. 

OBJECTIVES 

The  main  objective  of  this  project  is  to  examine  the  interaction  between  baroclinic,  mesoscale  eddies 
and  turbulence  using  a  large-eddy  simulation  (LES)  model.  Cases  will  focus  on  strong,  baroclinic 
waves  that  fonn  in  the  mixed  layer  along  surface  fronts  with  scales  of  a  few  km,  and  on  mesoscale 
eddies  that  are  imbedded  within  larger  scale  frontal  regions.  Our  goal  is  to  quantify,  understand,  and 
ultimately  parameterize  the  physical  processes  that  lead  to  lateral  mixing.  Simulations  will  help  guide 
field  experiments  planned  as  part  of  the  Lateral  Mixing  DRI,  and  provide  a  tool  for  understanding 
observations  in  the  analysis  phase  of  the  project. 

APPROACH 

High-resolution  simulations  of  baroclinic  instability  and  the  interaction  of  mesoscale  flow  with 
turbulent  mixing  are  conducted  and  analyzed  using  a  large-eddy  simulation  model.  Our  analysis 
centers  on  quantifying  and  understanding  the  mechanisms  by  which  small-scale  turbulent  structure 
develops  on  the  mesoscale  field,  the  physical  processes  and  balances  that  control  lateral  mixing  of  fluid 
properties  across  the  unstable  front,  and  the  transition  from  strongly  horizontal,  geostrophic  motion  on 
the  mesoscale  to  three-dimensional,  quasi-isotropic,  non-hydrostatic  motion  on  turbulent  scales. 
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WORK  COMPLETED 


Research  during  the  fouth  year  of  this  project  has  focused  on  completing  a  publication  focused  on 
simulations  of  baroclinic  instability  and  the  generation  of  turbulence.  Further  research  on  the  transfer 
of  energy  between  baroclinic  scales  and  turbulence  scales  was  also  conducted  and  is  reported  here. 

RESULTS 

A  key  question  in  our  analysis  of  small-scale  baroclinic  instability  concerns  the  fate  of  kinetic  energy 
generated  at  baroclinic  scales,  which  can  be  removed  from  the  system  by  either  a  gradual  cascade  of 
energy  through  progressively  smaller  scale,  quasi-two-dimensional  eddies  or  through  direct  production 
along  frontal  boundaries  of  three  dimensional  turbulence  that  decays  via  sub-grid-scale  viscous 
dissipation.  The  structure  of  kinetic  energy  spectra  as  a  function  of  horizontal  wave  number 

magnitude  k  =  ^ k2x  +  \cy  for  the  current  simulation  (from  a  two-dimensional  spectral  decomposition, 

averaged  over  wavenumber  vectors  of  approximately  equal  magnitude,  and  then  divided  by  k  to  obtain 
the  spectral  density)  also  suggests  that  energy  at  large  scales  is  not  removed  by  a  continuous  inertial 
cascade  process:  the  spectral  slope  is  about  -3,  more  indicative  of  an  energy-conserving  enstrophy- 
cascade  than  an  energy-cascade  regime,  and  there  is  enhanced  energy  and  a  nearly  flat  spectrum  at  the 
largest  wavenumbers,  between  scales  of  6  and  60  m  (Fig.  1).  A  rough  estimate  of  the  dissipation 
spectra  may  be  obtained  by  multiplying  the  kinetic  energy  spectra  (Fig.  1)  by  k~ ,  suggesting  a  modest 
peak  at  larger  wavenumbers  k  >  0. 1  m  ,  the  corresponding  variance-preserving  form,  obtained  by 
multiplying  by  a  third  factor  of  k,  would  indicate  that  the  dominant  contribution  to  the  integrated 
dissipation  is  from  wavenumbers  k  >  0.1  m  ,  or  at  scales  smaller  than  60  m  but  still  well-resolved  by 
the  3  m  grid  spacing. 


Wavelength  (m) 


Figure  1.  Two-dimensional  KE  spectra  calculated  as  described  in  the  text.  The  solid  line  is  from 
hour  84  and  dashed  line  is  from  hour  96.  Colored  lines  denote  various  spectral  slopes  as  labeled. 
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A  more  complete  description  of  the  spectral  scale  interactions  can  be  obtained  through  a  spectral 
decomposition  of  the  kinetic  energy  balance  (Frisch,  1995;  Capet  et  ah,  2008;  Molemaker  et  al.  2010), 
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where  ii(x,ky,zj  denotes  a  2-dimensional  spectrum  in  the  horizontal  direction,  asterisk  denotes 

complex  conjugate,  and  terms  are  defined  in  order  as  storage  tenn,  energy  flux  divergence,  pressure 
term,  buoyancy  tenn  and  dissipation  term.  For  this  analysis,  the  perturbation  velocity  and  the  terms 
that  enter  in  (1)  from  the  right-hand  side  of  the  momentum  balance  were  computed  and  integrated 
vertically  over  the  entire  model  domain,  and  separately  averaged  temporally  over  a  5  minute  time 
period,  with  the  product  then  providing  snapshots  of  the  spectral  energy  balance  at  two  different  model 
times,  hours  84  and  96,  representing  peak  wave  growth  and  reduced  growth  with  established  ring  or 
eddy  circulation  features,  respectively  (Fig.  2).  The  storage  tenn  in  (1)  was  calculated  by  computing 
the  difference  in  spectral  energy  over  the  5-minute  averaging  period  and  used  to  estimate  the  budget 
residual. 


(a)  (b) 


(c)  (d) 


Figure  2.  Domain  averaged  spectral  components  of  the  (a,b)  kinetic  energy  budget  and  (c,d) 
spectral  energy  flux  for  hours  (a,c)  84  and  (b,d)  96.  The  energy  flux  divergence  term  is  blue, 
buoyancy  term  red,  dissipation  term  green,  and  residual  is  black.  Spectral  components  in  (a,b)  are 

multiplied  by  wavenumber  k  to  preserve  variance. 
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The  integrated  spectral  fluxes, 


n(*)=  I  <j>(k)dk,  (2) 

k 

where  cp(k)  are  the  individual  terms  in  (1),  provide  a  measure  of  the  direction  of  energy  transfer  by 
each  term  in  the  budget.  For  example,  positive  values  for  the  integrated  energy  flux  divergence  term 
indicate  a  forward  transport  for  larger  wavenumbers,  with  energy  removed  by  dissipation  and 
numerical  effects  at  small  scales  (Fig.  2).  The  relatively  constant  values  for  the  dissipation  and 
residual  for  k  <  0.1  m'1  indicate  that  energy  loss  from  these  terms  occurs  only  at  the  small  scales  (k  > 
0.1  in1). 

Large-scale  baroclinic  waves  (k  <  0.01  m'1)  continue  to  gain  kinetic  energy  from  buoyancy  production 
generated  through  the  slumping  of  the  mixed  layer  front  through  hour  96  (Fig.  2a, b).  At  hour  84, 
large-scale  eddy  energy  is  lost  through  the  energy  flux  divergence  tenn,  even  though  the  overall  kinetic 
energy  is  still  increasing.  As  the  eddies  fonn  closed  circulations  at  hour  96,  the  energy  flux  divergence 
tenn  acts  to  redistribute  kinetic  energy  upscale  as  shown  by  the  negative  values  for  k  ~  0.003  m'1  and 
positive  values  for  k  ~  0.002  m'1.  Baroclinic  modes  clearly  dominate  the  energy  balance  at  large  scales 
while  turbulence  formation  and  dissipation  at  small  scales  close  the  budget.  There  is  almost  no 
indication  of  significant  energy  production  or  dissipation  at  wavenumbers  between  10'“  m'  and  10' 
m.  (Fig.  2a, b).  This  result  could  be  misinterpreted  as  an  indication  that  an  inertial  subrange  exists 
between  the  baroclinic  modes  and  turbulence  at  the  dissipation  scales.  However,  the  physical  structure 
of  the  flow,  with  energetic  turbulence  appearing  only  at  intensified  fronts,  and  the  reduced  strength  of 
eddies  at  intermediate  scales  indicated  by  the  steep  spectral  slope  (Figure  1)  suggest  that  turbulence  is 
gaining  significant  energy  primarily  through  non-local  spectral  transfer  rather  than  through  a 
continuous  forward  cascade. 

The  spectral  energy  budget  (8)  at  turbulence  scales  (Fig.  2a, b,  k>  0.1  m'1)  shows  a  peak  in  energy  flux 
divergence  associated  with  small-scale  instabilities  in  the  frontal  region.  At  these  scales,  perturbation 
kinetic  energy  is  lost  via  negative  buoyancy  production  when  vertical  mixing  destroys  stratification 
associated  with  the  sloping  fronts  and  the  mixed  layer  base.  Dissipation  is  a  maximum  near  the  grid 
scale  at  about  k  =  0.8  m'1.  However,  because  the  spectral  decomposition  (1)  of  the  energy  budget  is  not 
fully  compatible  with  the  numerical  discretization  and  because  the  model  grid  spacing  is  relatively 
coarse  for  simulating  turbulent  eddies,  there  is  also  a  fractional  residual  that  is  relatively  large  at  these 
small  scales,  where  the  budget  tenns  are  all  small.  The  residual  in  the  total  kinetic  energy  budget  is 
negligible.  The  flow-dependent  turbulence  closure  in  the  LES  model  allows  relatively  large  velocity 
fluctuations  at  small  scales,  with  a  steep  roll-off  of  the  kinetic  energy  spectrum  just  above  the  grid 
scale  (Fig.  1).  This  evidently  exaggerates  the  incompatibility  between  the  spectral  decomposition  (1) 
and  the  discretization,  leading  to  a  larger  residual  near  the  grid  scale  than  may  be  seen  in  hydrostatic 
circulation  models  that  spread  viscous  dissipation  over  a  wider  range  of  larger  scales,  such  as  the 
simulations  described  by  Capet  et  al.  (2008).  Note  in  this  context  that  a  related  incompatibility  is 
recognized  by  Capet  et  al.  (2008),  who  argue  that  their  near-zero  model  shear-production  tenn  can  be 
separated  ex  post  facto,  by  an  alternative  numerical  discretization,  into  a  positive  production  tenn  and 
a  balancing  negative  term  arising  from  the  model  discretization  that  is  presumed  to  represent  a  physical 
dissipation.  In  any  case,  our  main  conclusion  is  that  the  transfer  of  energy  from  the  larger,  baroclinic 
modes  to  small-scale  turbulence  is  a  more  direct  process  that  is  less  dependent  on  a  cascade  of  energy 
than  would  be  the  case  for  a  classical  shear  or  buoyancy  driven  turbulent  boundary  layer. 
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